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Abstract 

A useful connection between two-loop massive vacuum integrals and one-loop off- 
shell triangle diagrams with massless internal particles is established for arbitrary 
values of the space-time dimension n. 
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1. Many experiments testing the Standard Model and its extensions are sensitive to 
the values of two-loop corrections to physical quantities and require theoretical predictions 
for these contributions. Since evaluating massive two-loop diagrams is a tricky business, 
looking for non-trivial connections between different diagrams may be of certain interest. 

In a previous paper we noticed that, apart from some simple logarithmic terms, 
the finite part of the two-loop vacuum integral with three massive propagators involves a 
non-trivial function of the masses, which is exactly the same as a one-loop triangle with 
three massless propagators whose external momenta squared are equal to those masses 
squared. This suggests that there is a connection between these two seemingly unrelated 
diagrams. The purpose of this paper is to provide an explanation of this connection (for 
arbitrary values of the space-time dimension) and to explore some of its consequences. 

One of the important applications of the results for massive two-loop vacuum integrals 
is the small momentum expansion of two-loop diagrams with non-zero external momenta 
(see, e.g., in fj). Moreover, application of the general theory of asymptotic expansions 
H shows that analogous integrals also appear in the large momentum expansion and in 
the "zero-threshold" expansion (see in refs. @, |5|). Another application is the calculation 
of two- loop contributions to the p-parameter (see in || [7j). Furthermore, three- loop cal- 
culations in dimensional regularization M (with the space-time dimension n = 4 — 2s) 
require knowledge of the order s contribution to two-loop vacuum diagrams. As an ex- 
ample, we can mention recent calculations of three-loop corrections to the p-parameter 



and some other developments [10, 



The present paper is organized as follows. In section 2 we give definitions and differ- 
ent representations of the integrals considered. In section 3 we present some results for 
dimensionally regulated massless triangle diagrams. In section 4 we derive the connection 
between massive and massless integrals and use the corresponding results of section 3 to 
get the e-part of two-loop vacuum integrals with different masses. In section 5 we con- 
sider what the general results yield for the important case of equal masses. In section 6 
(conclusion) we discuss the main results of the paper. 

2. We shall discuss the following two types of Feynman integrals (see Fig. la,b): 

I(n; u 2 , , 3 K ml, m') = J J ^ _ {<f ^ _ ^ _ ^ , (1) 

f d n r 

J(n; , u ^\ P IvIpI) = j {{p2 _ rm(pi +rmr 2 K ( 2 ) 




Figure 1: Two- loop vacuum diagram (a) and one- loop triangle (b) 
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where p 3 = ~ (pi +^2)- To establish a connection between (HD and (||), we shall consider 
pf = m 2 , and we shall use dimensionless variables 

x = pl/pl = ml/ml , y = p\lp\ = m\jm\. (3) 

Below we shall omit the arguments m 2 and p 2 in the integrals / and J, respectively. 
Feynman parametric representations: 



/(n; 1/1,1/2, 1/3) = 1 vr — — / / / 



nr(i/j) 7 7 7 (a 1 a 2 + aiQ;3 + a2a3) ?1//2 (EajW 



2\Ei/i-n ' 



(4) 



w x -i-„ w /2 r(S Vi - n/2) /•/•/• n &<*i S(Z<Xi- 1) 

J(n;U U U 2 ,U 3 ) =1 71 I — — - 2 2 2 ^-»/ 2 > ( 5 ) 



II T(ui) J 7 7 (a^pj + aia 3 pl + a 2 a 3 p\f v -~ n l 2 ' 



Q;>0 



where J2 an d II denote the sum and the product from i — 1 to 3, respectively. Because of 
the delta functions, the integrations in (H) and (|5|) can be restricted to cti < 1. Neverthe- 
less, we prefer to write them in the present form, to simplify the discussion of parameter 
transformations we are going to describe below. 
Mellin-Barnes contour integral representations: 

= r ( , i /2) nrV.) W / / ds di x " v ' r( - s,rH) 

—loo —loo 

xr(n/2-u 1 -s)T{n/2-v 2 -t)r(v 1 + v 2 - n/2 + s + t)r (V f» - n + s + t) , (6) 



^n/2 jl-n / 2W2-£^ ^ lr f lc P 

r(n-E^i) nr(^i) (27ri) 2 / 7 

— loo —loo 

xT(n/2 - v 2 - v 3 - s)Y{n/2 -v x -v 3 - t)T{u 3 + s + t)T (J2 Vi~n/2 + s + t), (7) 

where the integration contours are chosen so as to separate the "right" and "left" series 
of poles of gamma functions in the integrand. By use of the residue theorem, the result 
for arbitrary n and z/j can be found in terms of hypergeometric functions of two variables 
(see 0,0, i). 

Uniqueness condition: When the i/j and n are related by v\ +V2+v 3 = n, a very simple 
result can be obtained from @ for such a "unique" triangle |13| (see also in |fL4|): 



J(n; ui,u 2 ,u 3 ) 



n/2 ;l-n 
71 1 1 



3 

n 



r(n/2 - Vi ) 



3. In the paper [15H, the following representation valid for arbitrary e was obtained: 



7T 



2 - £i i +2e r (i + e )r2(i_ e ) 1 r de r e ((v0' e - (*/0 



J(4 2 ' ;1,1,1) r(i-2e) J + + ® 
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We managed to generalize this result to the case when one of the z/s is arbitrary, as 



J(4-2e;l,l,l + 5) 



(p; 



T(l + 5) T(l-2e-8) 
1 rder £ ((2/0^-(^/0^) 



e + 5 7 (^2 + (! _ a- _ ^ + a-) 



1-e 



(10) 



At 5 = 0, eq. (fL0~l) gives (H), whilst for e = we get eq. (15) of |15| (see also in @). 
The expansion of J(4 — 2e; 1, 1, 1) in e is 

J(4 - 2e; 1, 1, 1) = tt 2 - i 1+2£ (^) _1 " £ r(l + e) {$ {1 \x, y) + e y) + (9(e 2 )} , (11) 

with 



$«(x,y) = - | 



d£ 



2/£ 2 + (1 - ac - 2/X + x 



ln- + 21nf) 

x J 



(12) 



vfr«(x,y) = - 



2/£ 2 +(l-£-?/)£+2 V x 



ln- + 21n£ 



In 



- ln(xy) 



(13) 

Due to the symmetry of the original triangle diagram with respect to all the external legs, 
the following combinations are totally symmetric in p 2 ,p|,p 2 : 



(p 2 )" 1 $ {1 \x,y) 



and 



(pi)- 1 (^(x^y) + | \n(xy) $«(x,y)) . (14) 



The integrals (p^)-([i3D can be evaluated in terms of polylogarithms, and the results 
can be found in []15| (eqs. (11) and (29), respectively). The results for $W were also pre- 
sented in fi~2| , U, |j~6| . Similar results for the triangle function in four dimensions (involving 
dilogarithms) were presented in |17|], while the e-part of massive triangle diagrams was 
considered in JTJ|. Here, we present the results for and \E , ^ 1 - ) in a different form: 



V) = 2^ ( 4Li2 (1 ~ Zl) + 4Lb (1 ~ Z2) + 4Li2 (1 ~ * a) 

+ In 2 Z\ + In 2 z 2 + In 2 z 3 + 2 In x In Zi + 2 In y In 22 1 , 

y) = --L {4Li 3 (l - ^r 1 ) + 4Li 3 (l - z 2 l ) + 4Li 3 (l - z^ 1 ) 
-4Li 3 (1 - zi) - 4Li 3 (1 - z 2 ) - 4Li 3 (1 - z 3 ) 
+41nxLi 2 (1 - z x ) + 41nyLi 2 (1 - z 2 ) 
— Yazi In z 2 In z 3 + In a; In zi ln(xzi) + lny In z 2 ln(yz 2 )j 

where the variables z^ and A are defined by 



(15) 



(X + x-y-lf 
Ay 



Z'2 



(X + y-l-x) 
Ax 



23 



(X + l-x-yf 
Axy 



(16) 



;i7) 
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X(x,y) = ^(1-x-y) 2 -Axy. 



The Zi are related to each other by ziz 2 z 3 = 1- The representations fllSQ and (jig) are 
equivalent to the ones presented in [15j in the region x + y < 1, as can be shown by 



transformations of the trilogarithms and dilogarithms [0. They are also valid in the 
regions x > y + 1 and y > x + 1. When A approaches zero, the Zi go to one, and hence 
all terms inside the braces vanishQ. The analytic continuations of (|15|) and (|T6|) into the 
region where x + y > 1 and x — 1 < y < x + 1, are obtained by adding some logarithmic 
terms to them, in such a way that the discontinuities along the branch cuts Zj < are 
cancelled. Using the formulae of fll9"l , the resulting expressions can be written as follows^: 



$«(x,y) 



v^A 2 " 
^ (1) (x,y) = In ( — J 



{Cl 2 (0i)+Cl 2 (# 2 )+ci 2 (0 3 )} 
2 



2 {Ls 2 (6 1 ) + Ls 2 (8 2 ) + Ls 2 (0 3 )} 



x,y) - 



Ls 3 (#i) + Ls 3 (# 2 ) + Ls 3 (63 



(19) 



7T 



irr < 20 > 



n o I + fl o / 1+35 — y \ /-l+x+y\ 

61 = 2 arccos [ — — ) , 8 2 = 2 arccos [ j= ) , 03 = 2 arccos [ — - — ) . (21) 



2^/y J ' ~ ~ '\ 2^ j 

Note that 9\ + 9 2 + 83 = 2-k. The log-sine integral is defined by 



2Jr\j 



Ls w (8) 



d8 In 



N-l 



8 

2 sin - 
2 



(22) 



and in particular, Ls 2 (8) = C1 2 (8). Eqs. (|T9|)-(pT|) are valid for all x, y such that A 2 <0. 

In the problem under consideration, the parabola defined by X(x,y) = is a special 
curve (see also in [pOf| ). It consists of the three segments \fx + ^fy = 1, ■sfx — yfy = 1 
and y/y — sfx = 1. On this curve, a simple result for arbitrary space-time dimension can 
be obtained using eq. @, 



J(4-2e; 1,1,1) 



7r 2 - £ i 1+2e T(e)r 2 (l - e) / x+y-1 y + 



lA=0 



+ 



X 



1+x-y 1 



2x 



(23) 



(pi)^ r(2-2e) ( 2.ny 2,/ 
which is clearly symmetric in p 2 , p 2 an d p 2 and can be written as 

- c rr r( rir (i 9~^ ) {(p^w + (papo^ 1 -* + ^(pd 1 -*} . (24) 

Pi V2 P3 1 l z — Ae ) 
Expanding ( p3|) in e gives 

2 In x 2 In y 



® {1) (x,y) 



A=0 



1 — x+y 1+x — y 



A=0 



ln 2 x— 41nx ln 2 y — 41ny 
1 — x + y 1+x — y 



(25) 



1 In Jig], the arguments of the L^'s and Lia's were chosen so as to vanish as x and y approach zero. 

2 Eq. (|19P corresponds to (4.15) of ffl). Representations similar to ( pQj ) were also considered in ^l) 
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4. Looking at the Mellin-Barnes representations ©— (0), it is possible to observe that 
r/. T/ „ „n ; i-n + v ?> ~ n / 2 ) T (^i + ^ ~ n/2)Y(u l + u 2 - n/2) 

l[n, Vi, V 2 , V 3 ) — 7T 1 — ; r— ; r— ; r 

r(z/i)r(z/ 2 )r(z/ 3 ) 

x J {TLvi - n; v 2 + v 3 - n/2, u x + u 3 — n/2, Vi + v 2 - n/2) . (26) 

The same relation can be obtained from the representations (§]) and (|5]) by first inverting 
and then rescaling the Feynman parameters ctj (i = 1, 2, 3) in (f|): 

-1 //— 1 1 //— 1 1 //— 1 

/ -7—/ // // //\ // -7—/ // // //\ t" 0^2 T Oi 3 fnn\ 

a i = —, a i= J r {a 1 ,a 2 ,a 3 )a i , F(a l ,a 2 ,a 3 ) = , (27) 

a • a" + a 2 + a 3 

where the scaling factor T has been chosen to restore the argument of the delta function 
to its original form. Due to the homogeneity of the integrand, the effect of the rescaling 
is to multiply it by a factor , which is precisely cancelled by the Jacobian associated 
with the change of variables a\ — > a" |22]. The result has the structure of (|5|). 

One's first impression may be that the relation (|26|) does not look very useful. In 
particular, in the case n = 4 — 2e, V\ — v 2 — v 3 — 1 it yields 

7(4 - 2e; 1, 1, 1) = vr 3 " 3 " i 1+2e T 3 (e) J(2 + 2e; e, e, e), (28) 

with some "strange" integral on the r.h.s. However, we can use the uniqueness relation 
(|8]) to transform the integrals J. Applying (^|) with respect to all three external legs of 
the three-point diagram gives: 

r(z/i)r(z/ 2 )r(i/ 3 )r(n - Ez/j) J (n; fx, v 2 , u 3 ) 

= r(n/2 - Vx)Y{n/2 - v 2 )T(n/2 - ^ 3 )r(E^ - n/2) 

x ^2 ) n/2-, 2 -,3( p 2 ) „/2-, 1 -,3( p 2 ) n/2-, 1 -, 2 J ( n; n/2 _ ^ n / 2 - y 2 ,n/2 - V 3 ) . (29) 

Combining (|26|) and (^9|), we get 



J(n; i*, i*, i*) = vr 3 "/ 2 - 2 - i 1 "" f n(m 2 ) n/2 ^) J (2Ei/ i -n; i*, i*, i*) . (30) 



\i=l 



Now, the powers V{ are the same on the l.h.s. and on the r.h.s. whilst the values of the 
space-time dimension are different. We shall call eq. (|30| ) a "magic" connection. It can 
also be derived by a change of variables in the Feynman parametric representation or in 
the Mellin-Barnes contour integrals. 

Let us consider what the "magic" connection gives for the most interesting case v\ = 
^2 = v 3 = 1 at different values of n. 

For n = 2 — 2e, we get 

1(2 - 2e; 1, 1, 1) = tT 3 * r 1+2e (m\m\m^ J(4 + 2e; 1, 1, 1). (31) 

So, the dimensionally-regularized triangle integral J(4 + 2e; 1, 1, 1) can be related to the 
two-dimensional integral 1(2 — 2e; 1, 1, 1) (also dimensionally-regularized, but note that 
signs of e are different!). Since both integrals are convergent as e — > 0, we can put e = 0: 

1(2; 1, 1, 1) = -i J(4; 1, 1, 1) = vr 2 r^ 2 $«(x, y). (32) 
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For n = 3 — 2e, we get 



7(3 -2e; 



7T 



3/2-3e :2+2e / 2 2 2\l/2-e 



[mim^m^j 



r(2e) 



r(3/2-e) 



J(3 + 2e; 1,1,1). (33) 



In this case, both integrals are three-dimensional, but (again!) the signs of e are different 
on the l.h.s. and on the r.h.s. Note that on the r.h.s. we also have a singular factor r(2e). 
So, to get the result for the singular and "constant" (in e) terms of 7(3 — 2e; 1, 1, 1), we 
need to know 7(3 + 2e; 1, 1, 1) up to the £-part. Using the representation @, it can be 
easily calculated, and we arrive at the following result: 

7(3 -2e; 1,1,1) = vr 4 " 2 " T 2 (l + e) j- + 2 - 4 In (m x + m 2 + m 3 )| + 0{e). (34) 

This corresponds to the result presented in [f23|, eq. (110). Note that simple results for 



three-dimensional triangles were presented in [^4 
For n = 4 — 2e, eq. (BUI) gives: 



7(4 - 2e; 1, 1, 1) = ^ i 1+2£ [m\m\m\)^ 7(2 + 2e; 1, 1, 1). (35) 

Using Feynman parameters, it is easy to show that 

7 (2 + 2e; 1, 1, 1) = -tt" 1 {7 (4 + 2e; 2, 1, l) + 7 (4 + 2e; 1, 2, l) + 7 (4 + 2 £ ; 1, 1, 2)} . (36) 

Now, we use the formula |12|] obtained by integration by parts |25|j , 

7(4 + 2s; 1, 1, 2) = {p\ ply 1 {-(pj+pl- p\) £ J(A + 2e; 1, 1, 1) 

+p 2 7(4 + 2s; 0, 2, 1) + p 2 7(4 + 2e; 2, 0, 1) - p 2 7(4 + 2e; 2, 1, 0)} (37) 

(note that the sign of e is different than in [Q). The integrals 7 with one of the v's equal 
to zero correspond to massless one-loop two-point functions. In such a way, we get 

{7 (4 + 2e; 2, 1, 1) + 7 (4 + 2s; 1, 2, 1) + 7 (4 + 2e; 1, 1, 2)} 
= (pIpIpI)' 1 {-eA(pl, pip 2 ,) 7(4 + 2e; 1, 1, 1) + p\{-p\ +p\+ p 2 )7(4 + 2e; 0, 2, 1) 

+pI(pI -Pl + Ps)7(4 + 2e; 2, 0, 1) + pl{p\ +p\- p 2 )7(4 + 2e; 2, 1, 0)} (38) 

where 

HpIpIpI) = 2p\p\ + 2p\p\ + 2p\p\ - {p\f - {p\f - {plf = -{plf A 2 (x, y) 



(39) 



is connected with the Kallen function. 

Thus, the connection can be written in the following symmetric form: 



7(4 -2e; 1,1,1)- , 

f 1 r(l + 2e) 

"^2+^1-25 r(i_ £ ) 



i vr 4 - 2e (m 2 m 2 m 2 )^ [(1 - e)(l - 25)]" 1 



x ■ 



7(4 + 2e; 1, 1, 1) A(m 2 , to 2 , to 



+r 2 (e) 



.TO 



2\e, 



— TO^-l-m^ + m^) + {m\) e {rrii — m^+m^) + (m^j^TO^-Hm^— to 



2\ei 



.(40) 



Using ([□]) to expand in e and keeping the terms up to the order e, we get 

7(4 - 2e; 1, 1,1) = | Tr^Cmi) 1 - 2 - T 2 (l + e) [(1 - e)(l - 2s)]- 1 

f 1 2 
x < — - (1 + x + y) H — (xlnx + ylny) — xln 2 x — yln 2 y + (1 — x — y) \nx\ny 

I e l e 

+^e(xln 3 x + yln 3 ?/) — |e(l — x — y) In a; lny (lnx + lny) 
- A 2 (a;,y)(l-e(lnx + ln?/)) $ (1) (a;, y) + e X 2 (x, y) ¥ 1] (x, y)\ + 0(e 2 ) . (41) 



The divergent and constant (in e) terms coincide with the result of Jl|, eq. (4.9). 

It could be noted that, using eqs. fl4D| ) and fl3"T|), we can write an exact relation between 
integrals I with different values of the space-time dimension: 

7(4 - 2e- 1, 1, 1) = [(1 - e)(l - 2s)]- 1 {-n' 2 ^ A(m 2 , m 2 , m 2 )/(2 - 2e; 1, 1, 1) 

+r 2 (e)(m 2 m 2 m 2 )~ e (mij^-mi+mj+m^ + tmj^mf-mj+mjj + ^j^mi+mj-mj) }. 

(42) 

Using (p0|) , an analogous relation can also be written for one-loop triangles, namely: 
J(2 + 2e; 1, 1, 1) = [plvlv^ e HpIpIpI) A* + ^; 1, 1, 1) 



T(l-e)T 2 (l + e) 1 



+2vr i+£ i 1 -* v r(1 + 2g) - [(PiPiMr + (?Wi)(pS)" + (paPa) W) e ] | -(43) 

This is a special case (iV=3) of an identity given in |26] relating one- loop iV-point integrals 
to (A/ - — l)-point ones. Note that at A = (A = 0) the first term in the braces on the 
r.h.s. disappears, and, changing e into (1 — e), we obtain nothing but the result 



5. In many realistic applications, the masses of internal particles are equal. For this 
case (jn% = m 2 = m 3 = m, x = y = 1), the integral representations (|1^) and ( Jl3| ) give 

1) = -2 J lnC, 1) = -2 / Ine In ( iZ f ± ^) ■ (44) 

We note that the first of these two integrals was also presented in f2Tfl . 

The angles ( pi]) are now all equal to 27r/3, so that eqs. (p0|) and (|l9f) are reduced to0 

$«(!,!) = -I Cl 2 (f ) = ^ Cl 2 (f ) , (45) 
*W(1,1) = ^ {-6Ls 3 (f)+4 ln3Cl 2 (f)-|7r 3 }. (46) 



Eqs. (|45|)-(f46|) can be expressed in terms of the generalized inverse tangent integral |fl9|] , 

1 f-l)^- 1 ) ln^V 

T ijv (z) = - (Li N iiz) - U N (-iz)) = z J d£— - A , (47) 



3 The representation in terms of Ck> (^) is well-known (see e.g. in || |2g|l). $W(1,1) can be also 
expressed as $W(1, 1) = | (|) - f tt 2 ] , see in ||. 
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whose Taylor series is 

Ti »W = ^g(^TTF = ^-3« + ^--. W<i- («) 

Using 

Ti 3 (73) = 73 £ y ( 2j { 1)3 = 0-570681635... (49) 

Cl 2 (|) = § Ti 2 (^) + i 7T In 3 = 1.014941606..., (50) 

L ^3 (?) = I Ti 3 (^) + f 1*3 Cl 2 (f) - i vr In 2 3 - ^ tt 3 = -2.144767213... (51) 
we find 

$ (1) (1, 1) = {l2 Ti 2 (^) + 7T ln3} = 2.343907239..., (52) 

1) = {-48 Ti 3 (^) + 7T In 2 3 + f 7T 3 } = 4.037576132... . (53) 



6. In this paper we have derived a useful relation (p0|) between two very different types 
of Feynman diagrams: two-loop massive vacuum integrals ([I]) and one-loop three-point 
functions (0). This "magic" connection is valid for any values of the space-time dimension 
and the powers of propagators. While the powers of propagators are the same, the massive 
and massless integrals related by (|30"D have different values of the space-time dimension. 
Nevertheless, using some additional transformations it is possible to relate dimensionally- 
regulated integrals considered around n = 4 (see eq. (|40"D ). However, we get different signs 
of e on the l.h.s. and on the r.h.s. (n = 4=p2e). As a result, some ultraviolet singularities 
of massive integrals can correspond to infrared singularities of massless triangles. This is 
not dangerous, since the magic connection relates diagrams of different type and does not 
introduce additional mixing of different singularities in diagrams of the same type. 

An important application considered in the paper is the result for the £-part of the 
two-loop massive vacuum diagram with different masses, eq. ([§!]). For the equal-mass 
case, a new transcendental constant is shown to appear. It is not excluded that it can be 
connected with an analytically-unknown constant occurring in three-loop calculations of 
the p-parameter (see also in p]|j ). 

Using integral representations fl9p-(|IUD, higher terms of the expansion in e (and in 5) 
can also be obtained. Moreover, the corresponding non-trivial functions will be the same 
in both cases, due to the magic connection. 
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